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Fractional Chern insulators are the proposed phases of matter mimicking the physics of fractional 
quantum Hall states on a lattice without an overall magnetic field. The notion of Floquet fractional 
Chern insulators refers to the potential possibilities to generate the underlying topological band- 
structure by means of Floquet engineering. In these schemes, a highly controllable and strongly 
interacting system is periodically driven by an external force at a frequency such that double tun¬ 
neling events during one forcing period become important and contribute to shaping the required 
effective energy bands. We show that in the described circumstances it is necessary to take into 
account also third order processes combining two tunneling events with interactions. Referring 
to the obtained contributions as micromotion-induced interactions, we find that those interactions 
tend to have a negative impact on the stability of of fractional Chern insulating phases and discuss 
implications for future experiments. 

PACS numbers: 73.43.-f, 05.30.-d, 67.85.-d, 71.10.Hf 


I. INTRODUCTION 

Ghern insulators^ form a primary and most widely 
studied class of more general topological insulators, 
proposed in condensed matter settings, in particular, 
heterostructures^^® and graphene.^’® They are character¬ 
ized by topological Bloch bands, that is, energy bands 
that give rise to a quantized Hall conductivity, when filled 
completely in a band-insulating state. The contribution 
of each band to the Hall conductivity is identified by 
an integer topological index, known as the Ghern num¬ 
ber. From this point of view topological bands can be 
thought of as generalizations of the Landau level.® How¬ 
ever, they can be realized in a broad variety of physi¬ 
cal settings independent of the requirement to have elec¬ 
trically charged particles coupled to an intense uniform 
magnetic field. The analogy between a Landau level and 
a Ghern band hints at the idea of fractional Ghern insu¬ 
lators (FGI).^®“^^ In this case, in addition to an isolated 
energy band characterized by a nonvanishing integral of 
the Berry curvature (which defines the Ghern number) 
one also needs strong particle interactions to form the col¬ 
lective states analogous to (and presumably richer than) 
the usual fractional quantum Hall states. 

It turned out that a powerful method to produce the 
desired topological bandstructures is Floquet engineer¬ 
ing. This form of quantum engineering is based on the 
fact that the dynamics of a time-periodically driven quan¬ 
tum system, a so-called Floquet system, is (apart from 
a periodic micromotion) captured by a time-independent 
effective Hamiltonian. Properties of the effective Hamil¬ 
tonian can be engineered by tailoring a suitable driv¬ 
ing protocol. Floquet engineering has been very suc¬ 


cessfully applied to quantum systems of ultracold atoms 
in periodically driven optical lattices.These sys¬ 
tems are particularly suitable for such control schemes 
due to their high degree of isolation from the environ¬ 
ment and versatile possibilities for controlling param¬ 
eters in a time-dependent fashion during the experi¬ 
ment. In recent years, several schemes have been pro¬ 
posed where periodic driving is employed for the realiza¬ 
tion of topologically nontrivial effective bandstructures 
in lattice systems, which in the absence of the driving 
are topologically trivial. These schemes can be divided 
into two classes: (i) methods working at high driving 
frequencies®®’^®’®^’®®’®®^®® and relying on averaging the 
driven Hamiltonian over a period, and (ii) methods work¬ 
ing at intermediate driving frequencies®®’®®^^® and going 
beyond the time-averaged description. For the latter, the 
term Floquet topologieal insulator has been coined. ^®’^^ 

The high-frequency schemes have been proposed and 
implemented in the context of ultracold atomic quan¬ 
tum gases in optical lattices®®^®®’®®’®^ and trapped ions.®® 
Here, static and time-periodic potentials are combined in 
such a way that the wavefunction acquires time-periodic 
relative phases on neighboring lattice sites. When av¬ 
eraged over one driving period, these phases resemble 
the nontrivial phases induced by a magnetic field. These 
schemes work at large frequencies, since it is assumed 
that the driving period T determining the phase mod¬ 
ulation is short compared to the tunneling time. The 
topologically nontrivial effective bandstructure of such 
a system has been recently probed in a square optical 
lattice.®^ 

The Floquet topological insulator schemes®®’^®’^® were 
originally proposed in the context of condensed matter 
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systems, considering irradiated graphene or semiconduc¬ 
tor heterostructures, and are based on a different princi¬ 
ple. Here, the effective Hamiltonian acquires new terms 
that describe tunneling between next-nearest neighbors 
and open a topological gap in the bandstructure. These 
new terms are related to second-order processes, where 
a particle tunnels twice during one driving period. Ob¬ 
viously they are significant only at intermediate driving 
frequencies, with the driving period being comparable to 
(or at most moderately shorter than) the tunneling time. 
This is the regime where particles undergo a significant 
periodic real-space micromotion in response to the driv¬ 
ing. Signatures of such Floquet topological bandstruc- 
tures have been observed in optical waveguides'^^ and 
with ultracold fermionic atoms in a shaken honeycomb 
lattice. 

As a natural next step, recently it has been proposed 
to stabilize a fractional topological insulator phase in 
such a Floquet topological bandstructure, referred to as 
the Floquet fractional Chern insulator.The scheme 
is based on fermions in a circularly driven honeycomb 
lattice, which acquires a topologically nontrivial effec¬ 
tive bandstructure.^® For their analysis, the authors em¬ 
ployed a high-frequency approximation^^^'^^ of the effec¬ 
tive Hamiltonian, including terms up to the second order. 
On this level of approximation, the nontrivial tunnel¬ 
ing terms opening the topological single-particle bandgap 
are included, but no corrections to the interaction terms 
appear. exact diagonalization of the approxi¬ 
mate effective Hamiltonian for small systems suggested 
that a topologically ordered state can be stabilized at a 
band filling of 1/3. 

The realization of such a topologically ordered many- 
body Floquet state is, however, challenged in various 
ways. One difficulty concerns the preparation of the state 
both in open condensed matter systems as well as in iso¬ 
lated cold-atom systems. Open Floquet systems assume 
steady states, which are generally quite different from 
equilibrium states'^®^®® and in an isolated system the de¬ 
sired state has to be reached by an adiabatic passage 
starting from the undriven ground states. A second con¬ 
cern is that excitations to higher-lying bands, which are 
not included in the tight-binding description, via multi- 
“photon” transitions may become relevant on the time 
scale of the experiment. Finally, also the impact of inter¬ 
actions beyond the second-order high-frequency approx¬ 
imation can be a relevant issue. Such “residual” interac¬ 
tions enter on two different levels, they cause heating and 
they lead to higher-order corrections to the approximate 
effective Hamiltonian. Interaction-induced heating cor¬ 
responds to processes that can be viewed as the resonant 
creation of collective excitations of the effective Hamil¬ 
tonian in high-frequency approximation. Such processes 
are not captured within the high-frequency expansion, 
they indicate that such an expansion cannot be expected 
to converge for an interacting system^^ and they are ex¬ 
pected to eventually drive the system towards an infinite- 
temperature regime.®^’®® A perturbative argument sug¬ 


gests that the rate of such detrimental heating will de¬ 
crease exponentially with increasing driving frequency 
(as the order in which the corresponding processes ap¬ 
pear increases with the driving frequency).®® Apart from 
heating, interactions will also lead to corrections appear¬ 
ing in higher orders of the high-frequency expansion, with 
the leading correction appearing in third order. 

In this paper we investigate the impact of leading in¬ 
teraction correction to the approximate effective Hamil¬ 
tonian. We do not address the issues of preparation and 
heating due to either multi-photon interband transitions 
or the resonant excitation of collective excitations. While 
it is rather clear that heating is a detrimental effect, it 
is an interesting question without an a priori obvious 
answer, whether the interaction corrections will tend to 
stabilize or to destabilize a Floquet fractional Chern in¬ 
sulator state. Even though the leading interaction cor¬ 
rections appear in third order only, they are still relevant 
for the required intermediate driving frequencies. First 
of all, if interactions are strong compared to tunneling, a 
third-order interaction correction can be comparable to a 
second-order kinetic term. And second, because the cor¬ 
rection has to be compared with the tiny manybody gap 
that protects the ground-state manifold of the effective 
Hamiltonian from excited states. The origin of the inter¬ 
action corrections to be investigated here is a significant 
real-space micromotion at intermediate driving frequen¬ 
cies. A particle at a certain lattice site will explore also 
neighboring sites during one driving period. This real- 
space micromotion generates new effective interactions 
at distances longer than those of the bare interactions 
characterizing the undriven model. 

For our study, we use exact diagonalization for small 
systems, taking into account the leading interaction cor¬ 
rection. In addition to a model of spin-polarized fermions 
with nearest-neighbor interactions (which includes the 
ground states found in Ref. 17) we also consider a system 
of spinless bosons with onsite interactions. The latter 
is particularly interesting for experiments with ultracold 
atoms in shaken optical lattices, and is typically governed 
by onsite interactions. In both models we find that taking 
into account interaction corrections tends to destabilize 
topologically ordered fractional Chern states. Thus the 
realization of Floquet fractional Chern insulator states 
seems rather challenging and it might be more promising 
to consider high-frequency schemes or schemes involving 
internal degrees of freedom.®^’®® 

The bulk of the presented material is split between two 
large sections, each subdivided into three subsections. 
Sec. H discusses the model, starting with the description 
of the lattice and the shaking protocol, and later proceed¬ 
ing to the Floquet analysis, the underlying single-particle 
problem and the role of interactions. Section HI focuses 
on results, and encompasses description of the numerical 
procedure, the obtained manybody bandstructures and 
quasihole spectra, characterizing the fractional Chern in¬ 
sulating states. Finally, we conclude with a brief sum¬ 
marizing Sec. IV. 
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II. MODEL 



FIG. 1. Honeycomb lattice as a triangular Bravais lattice with 
a two-site basis. The inequivalent sites A and B are marked 
with filled and open dots, respectively. Nearest neighbors are 
connected by the vectors <5^, while the set of next-nearest 
neighbor vectors Oj also defines the elementary translations. 

The model proposed by Grushin et al.^^ is based on a 
honeycomb lattice sketched in Fig. 1, and consisting of 
two intertwined triangular sublattices A and B. For the 
sake of further reference, Fig. 1 defines the vectors 


resented in the tight-binding form as a sum of the time- 
periodic kinetic and time-independent interaction parts, 

i?(t) =i?kinW+ffi„t. (4) 

The kinetic Hamiltonian reads 

3 

= — 'y ] y ] + h.c., (5) 

i^A fi—1 

here, the operator aj creates a spinless particle (fermion 
or boson) on lattice site i. The index /r labels the three 
distinct directions connecting any given site of sublattice 
A to its nearest neighbors belonging to sublattice B, and 
denoted by the shorthand label “i + /i” with i G A. Ex¬ 
tension to spinful particles, although not necessary for 
the purposes of the current presentation, is straightfor¬ 
ward. Equation (5) describes the nearest-neighbor tun¬ 
neling kinetics. The tunneling matrix elements along the 
three different directions acquire their time dependence 
due to the circular driving, and thus are time-periodic 
with uniformly distributed relative phases, i.e., 

Stt 

with ipf,= (6) 

Having in mind applications to ultracold atoms in a cir¬ 
cularly shaken optical lattice or graphene electrons irra¬ 
diated by circularly polarized light, we write 


6i=aey, 62\3 = T^ae^ + ^aey, ( 1 ) 

connecting a given site A to its three nearest neigh¬ 
bors (NN) on sublattice B, with a denoting the nearest- 
neighbor distance, and being the unit coordinate 
vectors. Likewise, the vectors ±a with 

Ctl — ^ 2|3 — zb ( 2 ) 

connect a given site of either type to its six next-nearest 
neighbors (NNN) belonging to the same sublattice. Note 
that NNN hopping transitions naturally separate into 
two classes. Transitions starting from a site A (resp. 
B) in the direction of aj correspond to counterclockwise 
(resp. clockwise) motion around the hexagonal cell, and 
the opposite statement applies to transitions along —a^. 
The vectors Oi and 02 are also taken to define the Bra¬ 
vais lattice r = niai -b n2a2, while = —ai — 02 is 
linearly dependent. 

The honeycomb lattice is subjected to a circular time- 
periodic force of fixed magnitude F, whose direction ro¬ 
tates in the {x-y) lattice plane at a constant frequency 
Lo, thus 

F{t) = F sin uit — F cos uit By. (3) 

This force results from either the circular shaking^^’^^’®^ 
of an optical lattice or, in the case of charged particles, 
from irradiation by circularly polarized light.A suitable 
gauge transformation^^’'^^ (corresponding to transition to 
the comoving frame of reference) restores the transla¬ 
tional invariance of the Hamiltonian, which is then rep¬ 


J{t) = J exp^ia sin uit), (7) 

with 


denoting the dimensionless shaking strength. 

In bosonic quantum gases contact interactions may be 
assumed, therefore, the interaction Hamiltonian is writ¬ 
ten in the standard Bose-Hubbard form 

Hint = y (9) 

ieA,B 

with the operator fii = ajfii measuring the particle num¬ 
ber on a given site. In contrast, for spinless fermions the 
onsite interaction term vanishes due to the Pauli exclu¬ 
sion principle, and repulsion between pairs of particles 
occupying neighboring sites must be taken into account, 
thus 

3 

HiH = ( 10 ) 

iGA ^1—1 

Therefore, depending on the particle statistics, we model 
the interactions using one of the two alternative forms 
given by either Eq. (9) or (10). Note also, that in both 
cases the interaction Hamiltonian depends only on den¬ 
sities, and thus remains static also in the presence of 
periodic driving. Our choice of treating spinless, that is 
spin-polarized, fermions is motivated not only by sim¬ 
plicity, but also by the observation of Ref. 17 that the 
relevant Floquet fractional Chern insulating states are 
ferromagnetic. 
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A. Floquet analysis 

As discussed in the Introduction, Floquet Chern insu¬ 
lators by construction require intermediate driving fre¬ 
quencies on the order of the tunneling strength or just 
moderately larger. This observation identifies the dimen¬ 
sionless inverse shaking frequency 


as the series-expansion parameter that classifies the suc¬ 
cessive contributions to the effective Hamiltonian, and 
in turn, to the ensuing physics. In further study we 
focus on the interval 0.1 ^ /3 ^ 0.5. At the lower 
limit of this range, one crosses over to the high-frequency 
regime where the effective Hamiltonian is adequately rep¬ 
resented by the time average of the driven Hamiltonian 
while the relevance of higher-order contributions fades 
away. In the opposite limit, for values of /3 exceeding 


one-half the representation of the effective Hamiltonian 
in terms of a /3-series expansion in no longer reliable. 

A powerful description of periodically driven quan¬ 
tum systems"^®^'*^’®^’®^ relies on the factorization of the 
quantum-mechanical evolution operator according to 

ti) = UF{t2) (12) 

thus separating micromotion described by the time- 
periodic unitary micromotion operator Upit) from the 
long-term dynamics captured by the effective Hamil¬ 
tonian Hp. A consistent high-frequency approxima¬ 
tion to the effective Hamiltonian Hp is constructed in 
Refs. 45, 47, and 62 following different approaches and 
resulting in equivalent®^ series representation in powers 
of w-i 

Hp = H^j^'>+ hP+hP+ ■■■ , (13) 

where 



m—1 


M3) ^ 1 ^ 1 [H-m', [Hm'-m,Hm]] - , [H-m'-n., Hm]] 

^ 2{fiujY _ rrfl ‘iitnoY , mm' 

1 m,m'— 1 

m'^m 


(14a) 

(14b) 

(14c) 


are expressed in terms of Fourier components of the 
driven Hamiltonian 

OO 

H{t)= Y. (15) 

m— — oo 

It is useful to note the hermiticity condition H-m = Hj^. 

For the purposes of the current application, the general 
expressions (14) can be simplified in two aspects. First, 
we note that the time dependence of the driven Hamil¬ 
tonian stems entirely from the time-dependent hopping 
amplitudes described by the Fourier expansion 

OO 

Mt)= E JJm(a)e-™^-e™‘"*, (16) 

m— — oo 

with Bessel functions of the first kind (and order to), 
denoted here by Anticipating moderate shaking 

strengths a [see Table I] we note that the Fourier series 
of the Hamiltonian converge very fast since contributions 
originating from the mth harmonic with to > 0 include 
a prefactor ^/[^(a) ^ (see Fig. 2). For this rea¬ 

son, we truncate the Fourier expansion and include only 
the terms with |to| ^ 1. Consequently, non-zero Fourier 


components of the driven Hamiltonian are 


Ho = Hint - EE JjJo{ct) -l-fij, (17a) 

iGA M 

H±i = T EE JJi(a) 

iGA M 

(17b) 


aj+^ai - a\ai+^ 


where the interaction Hamiltonian Hint is given by either 
Eq. (9) or Eq. (10) for bosonic and fermionic systems, 
respectively. 

Secondly, focusing on the interplay of micromotion 
and interactions we take into account the leading third- 
order interaction correction while neglecting purely ki¬ 
netic terms of third order. Thus, the expansion of the 
effective Hamiltonian considered in our work reads 


H^P = Ho, 

_ 

^ 2(fia;)2 


Hi, Hi 


Hi, 


Hint, Hi 


h.c 


(18a) 

(18b) 

(18c) 






















5 



FIG. 2. Comparison of the importance of various Fourier 
harmonics. The shaded background shows the range of the 
dimensionless shaking amplitudes a considered in the calcula¬ 
tions. The red lines depict the behavior of J7 o(Q!) and 
setting the scale of, respectively, the static components of the 
Hamiltonian and the first harmonic. The dashed blue lines 
show the behavior of the (omitted) subleading contributions 


B. Single-particle spectrum 

The first-order contribution to the effective Hamilto¬ 
nian given by Eq. (18a) is obtained by averaging the 
driven Hamiltonian F?(t) over a period. This leads to the 
result given by Eq. (17a) and featuring the well known 
renormalization of NN hopping amplitude according to 
the prescription J —>■ Jl7o(ct), which has been extensively 
exploited to demonstrate a number of interesting physi¬ 
cal effects in cold atom settings, including the superfluid- 
Mott insulator phase transition. 

The second-order Floquet contribution is expressed by 
the commutator in Eq. (18b), and physically corresponds 
to two subsequent tunneling processes of a particle dur¬ 
ing one driving period. The basic commutation relation 
[dlae^al^an] = 5(,md\an - 5knd\ndt is valid for particles 
of either statistics and combines two successive NN tun¬ 
neling events into an effective NNN transition. Thus, 

Hf = - E (19) 

«b» 

with the sum running over all next-nearest neighbor con¬ 
nections, and 

=±iv^Wf(a). (20) 

The plus (minus) sign applies to clockwise (counterclock¬ 
wise) transitions around the hexagonal unit cell. Note 
that NNN transitions in Eq. (19) are characterized by 
purely imaginary amplitudes. 

The role of the second-order contribution to the ef¬ 
fective Hamiltonian in Floquet engineering is to open 
topological gaps in the single-particle spectrum. This 
produces the Haldane model^ using driven honeycomb 
lattices.Since the possibility to stabilize the 
FCI phases is the ultimate question of this work, one 


is interested in starting from a favorable single-particle 
bandstructure characterized by relatively flat bands (at 
least in a part of the Brillouin zone) that are separated 
by large energy gaps. As noted previously,this require¬ 
ment leads to a constraint relating the shaking amplitude 
a to the inverse frequency /3. The detailed discussion of 
this matter is delegated to the Appendix A, with the 
conclusion that the relative strengths of NNN and NN 
transition amplitudes must be close to the ratio 


Jo(a) 4-\/6 


( 21 ) 


Thus, for any given value of the inverse shaking strength 
/3, the constraint (21) defines the corresponding “nom¬ 
inal” value of the driving strength oq which optimizes 
the single-particle bandstructure. Let us stress, how¬ 
ever, that thus defined value ao should be regarded more 
as representative rather than precisely defined optimal 
value. We verified that moderate deviations of the value 
of a from ao do not change the results qualitatively. 
Thus, focusing on a single value helps to reduce the di¬ 
mensionality of the parameter phase space. 


C. Role of interactions 


A cornerstone of the present contribution is the argu¬ 
ment that the real-space micromotion couples to the par¬ 
ticle interactions through Eq. (18c), which leads to the 
generation of new (and modification of existing®®) inter¬ 
action terms, and in particular, influences the formation 
and stability of FCI phases. Possible interplay of mi¬ 
cromotion and interactions was previously discussed®®^®® 
basing on approximation schemes related to the Magnus 
expansion.®® The obtained terms are proportional to the 
inverse driving frequency and typically offer a clear phys¬ 
ical interpretation (such as the density-assisted tunnel¬ 
ing). However, as discussed in Ref. 47 as well as Ref. 45, 
these terms do not influence the spectrum within the or¬ 
der of the approximation. 

In the present context of circularly shaken honey¬ 
comb lattices, micromotion-induced interaction correc¬ 
tions were analyzed and the physical nature of the ad¬ 
ditional terms was identified in Ref. 47. Following this 
work and specializing to moderate driving amplitudes, 
such that Fourier components \m\ ^ 1 are sufficient, we 
observe that the overall strength of these terms is set by 
the prefactor 


UJ^J^ia) 

{hwY 


UP^Jl{a) = pU. 


( 22 ) 


Here we introduced the dimensionless quantity rj = 
0^Ji{a) as a natural measure of the relative strength 
of micromotion-induced interactions with respect to the 
bare onsite repulsion. The five additional contributions 
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for bosons with contact interactions read: 

^ U^n, {ui - 1) (23a) 

i 

+ 477 [/ fiihj (23b) 

(ij) 

+ 2r]U'^alalajaj (23c) 

(b> 

- ^rjU al{Anj - fii - nk)ak (23d) 

{ijk) 

+ 577C/ ^ (a]a]aia/c + h.c.^ . ( 23 e) 

{ijk) 


Here, 2 ; = 3 is the coordination number (the number of 
nearest neighbors), and the sums are taken, respectively, 
over the lattice sites i, all directed NN links (ij), and all 
directed three-site strings {ijk) with i and k being next- 
nearest neighbors connected via an intermediate site j. 
The physical interpretation of the obtained terms is as 
follows: (a) reduction of the onsite interaction strength, 
(b) nearest-neighbor density-density interaction, (c) pair 
tunneling, (d) density-assisted tunneling between NNN 
sites, and (e) cotunneling of pairs of particles into (from) 
a given site from (into) two distinct nearest neighbors. 
While for fermionic particles the corresponding analysis 
becomes more involved and is not presented here, we em¬ 
phasize that the general comparison of the relative impor¬ 
tance of micromotion-induced interactions and the defi¬ 
nition of rj remain unchanged except for the replacement 
of onsite repulsion energy U with NN repulsion energy 

In general, the overall scale of the micromotion- 
induced interactions (i. e., the value of 77 ) is not an in¬ 
dependent parameter but is set by the driving frequency 
(/3) and strength (a). Since the preferred value of the 
driving strength oq is in turn fixed by the flat-band con¬ 
dition, all in all, we have to consider a one-dimensional 
cut through the apparently three-dimensional parameter 
space (a, p, rf). The characteristic range of values is sum¬ 
marized in Table I. 


TABLE I. Numerical values of parameters used in calculation. 


p = 0.1 

ao = 1.814 

7 = 0.003 

p = 0.2 

ao = 1.391 

0 

d 

II 

P = 0.3 

ao = 1.153 

7 = 0.021 

d 

II 

ao = 1.004 

7 = 0.031 

10 

d 

II 

ao = 0.900 

7 = 0.041 


As explained previously, the eligible values of the dimen¬ 
sionless inverse shaking frequencies P are constrained to 
an interval consistent with the physical setting. One can 


see, that the corresponding values of the shaking strength 
ao are also restricted to an interval where: (i) the single¬ 
particle bands do not collapse as the zeros of the zeroth- 
order Bessel function are avoided, and (ii) the employed 
truncation of the Fourier series is indeed valid. 

Turning to the last column displaying the values 
of 77 , one observes that the relative contribution of 
micromotion-induced corrections to the overall strength 
of interactions is limited to a few percent. However, 
these low values should not be taken as a proof that 
micromotion-induced interactions may be neglected. In 
fact, in the context of the stabilization of FCI phases, 
their importance should be measured with respect to 
manybody gaps separating the ground-state manifold, 
and judged by their impact on the FCI stability diagram. 
As we will see shortly, micromotion-induced interactions 
may indeed be significant and tend to destabilize frac¬ 
tional Floquet states. 


III. RESULTS 
A. Numerical procedure 

To gain specific insights into the impact of micromo¬ 
tion on the formation and stability of FCI phases we per¬ 
formed a numerical exact-diagonalization study of peri¬ 
odically driven honeycomb lattices. The most stable frac¬ 
tional quantum Hall-like states are expected to form at 
Ailing fractions k' = ^ for fermions and k' = ^ for bosons. 
Therefore, our numerical modeling mostly focused on fi¬ 
nite systems of Np = 8 spinless fermions (resp. bosons) 
moving on a lattice containing Ni x N 2 = 6 x 4 (resp. 
iVi X A 2 = 4 X 4) elementary two-site cells. While com¬ 
parisons with results obtained for larger systems were 
used for consistency checks and validation of the general 
conclusions, the bulk of calculations was performed on 
moderate eight-particle systems. Keeping the system size 
relatively small allowed for looping over a dense grid in 
the parameter phase space defined by the aforementioned 
inverse driving frequency P and the interaction strength 
U or V, as well as the auxiliary phases 71 and 72 intro¬ 
duced through twisted boundary conditions^® in the two 
primary lattice directions. Twisting the boundary condi¬ 
tions in the ith direction by a phase factor e‘'''’ represents 
the insertion of the dimensionless flux 7 ^, and leads to the 
variation of the calculated energy levels, which is com¬ 
monly referred to as the spectral flow. The formation of 
the topological order is signaled (however, not rigorously 
proven) by the formation of the ground-state manifold 
(GSM) consisting of quasidegenerate states that re¬ 
arrange under the spectral flow but remain isolated from 
the other states by a finite manybody gap.^^^s 

In view of the translational invariance, restored in the 
finite system by the cyclic boundary conditions, our cal¬ 
culations are performed in the reciprocal (or quasimo¬ 
mentum) fc-space. The finite geometry of the lattice im¬ 
poses a discreet grid of Ni x N2 permissible values in 
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a Brillouin zone (BZ) for both single-particle and total 
quasimomenta. The diagonalizations can be performed 
separately at each total quasimomentum k thus signifi¬ 
cantly reducing the size of the problem. To quote a spe¬ 
cific example, in the 8-fermion system the dimensionality 
of the Hilbert space reduces from 24!/(16! 8!) ~ 7.35 x 10® 
to just above 3 x 10^ when translational invariance is 
taken into account. The fc-space points belonging to a 
single Brillouin zone are labelled by integer pairs (^ 1 ,^ 2 ) 
with ki e [0, Ni — 1] and ^2 G [0, N 2 — 1] or, equivalently, 
by a single index K = ki + Nik 2 with K G [0, fVi7V2 — 1]. 
The total quasimomentum sectors at which the FCI 
states will form are predicted by a simple counting rule^^ 
whose validity in the present case is indeed confirmed by 
the actual numerical results. Thus, for eight fermions, 
the FCI states form at K = {0,2,4}, while for eight 
bosons one finds two states at TF = 0. 

Following previous works that demonstrated the 
potential existence of FCI phases in various lattice 
models^'^^^^ with strong particle interactions, we use the 
customary band projection technique (see, e.g., reviews 
in Refs. 15 and 16 for an extensive discussion). Thereby, 
only processes in the relevant single-particle band are in¬ 
cluded, while the interband scattering events as well as 
processes related to the remaining band are neglected. 
Note however, that we do not flatten the bandstructure 
which is another standard tool of the trade.The rea¬ 
son for this choice of techniques is the following: In con¬ 
trast to models based on energy bands featuring weak 
dispersion throughout the Brillouin zone, we deal here 
with single-particle bands that are in general dispersive, 
however, the upper (resp. lower) single-particle band de¬ 
velops a relatively fiat section at the bottom (resp. top) 
where also the Berry curvature is the largest. In this sit¬ 
uation, band flattening would lead to an essential distor¬ 
tion of the model and therefore must be avoided. When 
modeling fermionic systems we adopt the approach of 
Ref. 17 and consider a gas at the density of 4 particles 
per 3 elementary cells. This produces a completely filled, 
and thus inert, lower band plus a partially filled upper 
band with particles predominantly concentrating in the 
fiat section. Therefore, projecting onto the upper band 
we are able to focus on the physics at the effective filling 
factor ly = ^. Proceeding to bosonic systems we envi¬ 
sion a gas at the density of one particle per 2 elementary 
cells, and employ projection onto the lower single-particle 
band. 


B. Manybody topological gap 

As discussed previously, the parameter phase space 
to be explored is essentially two-dimensional, and is 
spanned by the inverse driving frequency j3 = J/hio^ and 
the particle repulsion strength. It is natural to measure 
the interaction energies with respect to the basic hop¬ 
ping amplitude, thus introducing the corresponding di¬ 


mensionless quantities 

U V 

u = —, and V = —, (24) 

for bosonic and fermionic systems, respectively. 


micromotion omitted micromotion included 



FIG. 3. Typical energy spectra of the driven honeycomb 
lattice with micromotion-induced interactions ignored [panel 
(a)], and taken into account [panel (b)[ plotted versus the aux¬ 
iliary flux 71 inserted along the first (commensurate) lattice 
direction. Here, the energies of manybody states in the eight- 
fermion system numerically calculated at /3 = 0.3, v = 5, and 
72 = 0 are shown. Colored (resp. black) lines denote the states 
in the ground-state manifold (resp. other) total quasimomen¬ 
tum sectors, and two lowest-energy states of each quasimo¬ 
mentum sector are included. The dark red arrow illustrates 
the definition of the local manybody gap referred to in the 
text. 

Figure 3 shows the calculated energy spectra of eight- 
fermion manybody states in the form of spectral flows 
along the first (commensurate) direction. The three 
states belonging to the expected FCI total quasimomen¬ 
tum sectors K = { 0 , 2 , 4 } are plotted in color while the 
remaining states are plotted in black. The insertion of 
the artificial flux 71 defining the twisted boundary con¬ 
ditions makes the three states in the GSM interchange 
and reconnect to their partners at the opposite bound¬ 
ary of the manybody Brillouin zone (MBZ). The plot is 
obtained by setting /3 = 0 . 3 , v = 5, and 72 = 0 . In 
the left panel, the micromotion strength 77 is artificially 
set to zero (as would happen in a calculation omitting 
the presence of real-space micromotion). The nearly de¬ 
generate states in the GSM stay protected from the ex¬ 
cited states by a manybody gap, and their total Chern 
number obtained from sampling over the whole MBZ 
(71 172) G [ 0 , 27 r) X [ 0 , 27 r) sums up to unity with numerical 
precision. The right panel shows the realistic situation 
where coupling between micromotion and interactions is 
duly taken into account up to the third order. While the 
topological nature of the lowest-energy manifold and the 
manybody gap persist, its width is reduced and the en¬ 
ergy spread of the states in the GSM is now considerably 


















larger. These observations serve as an early hint that the 
interplay of micromotion and interactions is significant, 
and may indeed have a detrimental role on the stability 
of the Floquet fractional Chern insulating phases. 

Referring to Fig. 3, we also take the opportunity to give 
a precise definition to the notion of the manybody gap to 
be used in presentation of further results. Since calcula¬ 
tions are performed looping over the manybody Brillouin 
zone we first define the local (that is, corresponding to 
fixed values of 71 , 2 ) dimensionless manybody gap 

= 7 • 

(25) 

as the energy difference between the highest state belong¬ 
ing to GSM and the lowest-lying state outside the GSM 
measured in units of J. Thus, A( 7 i, 72 ) is positive when 
an isolated ground state manifold is formed and negative 
when the states in the ground-state manifold mix with 
the remaining states. The overall manybody gap is then 
obtained by minimizing the local manybody gap over the 
entire MBZ, thus 


A = min A( 7 i, 72 ). 

MBZ 


(26) 



FIG. 4. The phase diagram showing the manybody gap 
A as a function of the inverse shaking frequency /3 and the 
dimensionless interaction strength u for eight-boson system. 
Micromotion-induced interactions are omitted in panel (a) 
and taken into account to the third order in panel (b). The 
full black lines delimit the regions of positive manybody gaps 
(shown in red shades) and differ considerably between the 
panels. The white dashed line indicates the parameter regime 
where the interaction strength U is equal to the single-particle 
bandgap. Well below this line mixing between single-particle 
bands becomes small. 


The complete phase diagrams in the j3-u (resp. f3-v) 
plane for bosons (resp. fermions) occupying the lower 
(resp. upper) single-particle band at 1/3 (resp. 1/2) fill¬ 
ing are shown in Figs. 4 and 5. In both instances we map 
out the behavior of the manybody gap as a function of 
the governing parameters, /3 and u (resp. v). The red 
areas correspond to the presence of a positive manybody 
gap implying that the states comprising the GSM are 



FIG. 5. The same as Fig. 4 but for eight-fermion system with 
dimensionless NN repulsion energy v. Note that the scale of 
the colorbar, and correspondingly the typical manybody gaps 
are significantly smaller than for bosons. 


isolated from the remaining ones everywhere in the MBZ 
and thus the two manifolds do not mix. On the contrary, 
blue areas correspond to a negative manybody gap. In 
this situation, at least one state from the GSM attains 
higher energy value than the lowest state not belonging 
to the GSM, and the two manifolds overlap somewhere in 
the MBZ. The full black line separates phase space areas 
characterized by manybody gaps of different signs and 
thus serves as a first criterion for the possibility of forma¬ 
tion of the Floquet fractional Ghern insulator. Of course, 
this criterion takes into account only the aspects of the 
energy spectra and must be further supported by, e. g., 
information obtained from excitation spectra discussed 
in the following subsection. Nevertheless, already these 
phase diagrams indicate that inclusion of micromotion- 
induced interactions typically pushes the phase boundary 
upwards in the phase diagram. That is, stronger inter¬ 
actions will be needed to stabilize FGI phases at equal 
other conditions. 

Bosonic lattice systems seem to be much more promis¬ 
ing candidates for the stabilization of FGI phases than 
fermionic ones. Indeed, they find a realization in quan¬ 
tum gas experiments where strong and tunable onsite in¬ 
teractions are possible, and Fig. 4 reveals that in bosonic 
systems large manybody gaps on the order of the tun¬ 
neling amplitude J can be obtained. In contrast, for 
fermionic systems nearest-neighbor repulsions are harder 
to tune into strongly-interacting limit (in graphene the 
strength of NN repulsion is estimated^^ to be around 2J), 
and the overall scale of attainable manybody gaps is also 
smaller, as demonstrated by Fig. 5. 

Viewing FGI states as direct analogues of Laughlin’s 
fractional quantum Hall states, one would ideally prefer 
working in the regime where the fractional states are in¬ 
duced by strong interactions mixing single-particle states 
within a single topological (Ghern or Landau) band and 
not among several bands. (Recent results suggest, how¬ 
ever, that weak interband mixing is not detrimental for 
the formation of fermionic v = ^ FGI state^^’^^). This 
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corresponds to the situation where the single-particle 
bandgap exceeds the characteristic interaction strength. 
In Figs. 4 and 5, the white dashed lines indicate the 
regime where the corresponding interaction parameter, U 
or V, is equal to the single-particle bandgap; well below 
these lines mixing between single-particle bands is small. 
Obviously, the approximation based on the band projec¬ 
tion is not justified in the considered regime. Moreover, 
the interaction strengths required to stabilize FCI phases 
in the current model also generally exceed the above- 
stated limit. These observations suggest that it must be 
quite challenging to explore FCI phases within the con¬ 
sidered model. A window of opportunity might persist 
at the higher end of the considered range of driving fre¬ 
quencies (equivalently, at the lower end of the range of 
/3’s). Although the approximation focusing on just one 
single-particle band is not fully justified here, the effects 
of band mixing were reported not to be critical.In 
this limit, notably weaker interaction strengths are re¬ 
quired to stabilize FCI phases [see Figs. 4 and 5], and 
moreover, problems caused by heating due to the reso¬ 
nant excitation of collective excitations become smaller. 


C. Quasihole spectra 


micromotion omitted micromotion included 



K K 


FIG. 6. Quasihole excitation spectra for eight-boson system 
calculated at the inverse dimensionless driving strength P = 
0.1. Formation of an isolated manifold of 4x 16 states is visible 
for strong interactions with u = 5, and gradually closes when 
interaction become weaker. 

Finally, let us also take a brief look at the quasihole 
excitation spectra, corresponding to the removal of a par¬ 
ticle from the ground state.Here, we are driven by 
twofold motivation. On the one hand, these spectra serve 
as an additional identification check confirming that the 
obtained phases display behavior characteristic of frac¬ 
tional Chern insulator. Thus, one expects to observe the 
formation of an isolated low-energy manifold of quasihole 


states whose number is predicted by counting rules based 
on the generalized Pauli principle.On the other hand, 
one is interested in the the evolution of the size of the 
gap separating the low-energy manifold from the remain¬ 
ing states as a function of the computational parameters. 
In this way, we also obtain complementary information 
about the interaction strengths that are sufficient to in¬ 
duce the FCI phases. 

Qualitatively similar results are obtained for both 
fermionic and bosonic systems. Focusing on the lat¬ 
ter, we perform exact diagonalizations for systems of 7 
bosonic particles moving on a lattice consisting of 4 x 4 
elementary cells. This is one particle less than would 
correspond to the exact 1/2 filling of the band, and can 
be interpreted as an introduction of a hole into the pre¬ 
viously studied system. In this specific case, the (1,2)- 
admissible counting rule^^ predicts the formation of an 
isolated manifold with 4 states per each total quasimo¬ 
mentum sector, implying the total of 64 states. The ac¬ 
tual numerical results obtained for /3 = 0.1 and varying 
values of the interaction strength are shown in Fig. 6. In 
all panels, the abscissa axes enumerate the total quasi¬ 
momentum sectors indexed by the integer values K run¬ 
ning from zero to 15 = 4x4—1, and the ordinate axes 
display the energies of the lowest manybody states ob¬ 
tained from the exact-diagonalization calculations. The 
left (resp. right) column of three plots corresponds to 
omitted (resp. included to the third order) micromotion- 
induced corrections to the effective Hamiltonian. The 
strength of particle interactions u = U/J is decreasing 
from the top to the bottom. At the strongest considered 
interactions, u = 5 (top row) one sees a clear gap sepa¬ 
rating the lower group of exactly four states per quasi¬ 
momentum sector from the rest of the spectrum. For 
weaker interactions (u = 3 in the middle row) the gap is 
barely discernible in the absence of micromotion-induced 
interactions and is obliterated when micromotion is taken 
into account. Finally, the in bottom row corresponding 
to weak interactions with u = 1 there is no visible gap 
indicating the absence of FCI phase. In general, the in¬ 
formation provided by the quasihole spectra is broadly 
compatible with that obtained from the ordinary many¬ 
body spectra, however, the constraints placed on interac¬ 
tions being “sufficiently strong” are even more stringent, 
thus further contributing to the pessimistic outlook on 
the feasibility of fractional Chern insulating states in the 
studied system. 


IV. CONCLUSIONS 

To summarize, we address the issue of the stabilization 
of Floquet fractional Chern insulator states for strongly 
interacting particles (bosons or fermions) on a time- 
periodically driven honeycomb lattice. In this system, 
the necessary topological single particle bands are formed 
due to the next-nearest neighbor transitions, which is 
a second-order effect corresponding to two consecutive 
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tunneling events during a single driving period. This 
requires sufficiently low driving frequencies and necessi¬ 
tates the consideration of further expansion terms beyond 
the averaging of the Hamiltonian used in high-frequency 
schemes. The third-order terms describe the coupling of 
real-space micromotion and interactions. In strongly in¬ 
teracting systems, the importance of these terms is, in 
general, comparable to that of the second-order contri¬ 
butions. The prefactor features not only the expansion 
parameter but also the onsite repulsion energy U 

for bosonic systems or, alternatively, nearest-site repul¬ 
sion energy V for fermionic systems. Within simulations 
of small systems, the coupling of micromotion and inter¬ 
actions turns out to be both significant and detrimental 
to the formation of quantum-Hall like states. Thus, the 
realization of Floquet fractional Chern insulator states 
seems rather challenging. 


[0,27r) X [0,27r) covering the rhombic Brillouin zone. In 
the case considered in the main text, the NNN hop¬ 
ping amplitude is purely imaginary, corresponding to 
—)■ 7r/2, and therefore 


H{k) 


-J 2 f{k) -Jig*{k) 
-Jigik) J 2 f{k) 


with the unchanged definition for g{k), and 


f{k) = 2 [sin ki + sin ^2 — sin(A:i -I- ^ 2 )] • 


(A4) 


(AS) 


While the resulting bands are not globally flat, they fea¬ 
ture flat sections that cover BZ regions characterized by 
large Berry curvatures. This fact forms the basis for 
the anticipation^^ that the Haldane model might support 
FCI phases for strongly interacting particles. 
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Appendix A: Haldane model 


The Haldane modeb’^®’^^’^®’^^^^® sets a paradigmatic 
example of the topological bandstructure supported by a 
simple two-band setup. This model features a honeycomb 
lattice (see Fig. I) with nearest-neighbor hopping de¬ 
scribed by a real amplitude — Ji plus next-nearest neigh¬ 
bor hopping described by a complex amplitude — J 2 e“^ 
(resp. —J 2 e““^) in the counter-clockwise (resp. clock¬ 
wise) direction. In the quasimomentum representation, 
the Hamiltonian matrix reads 


Hik) = 


-^2/(fc,^) 

-Jig{k) 


-Jig*{k) 

-J2f{k,-if>) 


(Al) 


with 

3 

g{k) = ^ (A 2 ) 

3 

f{k, (^) = 2 ^ cos (fc • aj - If) 
i=i 

= 2 [ COs(fci — (fi) + cos(fc2 — f) 

+ cos(/ci + k 2 + </?)] • (A3) 



FIG. 7. Dispersion of the upper band of the Haldane model 
drawn along the diagonal of the rhombic BZ. 


The formation of the flat segments in the energy 
bands is best visualized focusing on a one-dimensional 
cut through the rhombic BZ along its diagonal (obtained 
by setting ki = ^ 2 ). This line connects the zone center 
at (0,0) to the equivalent point at (27r, 27r), and passes 
through both inequivalent Dirac points at (27r/3,27r/3) 
and ( 47 r/ 3 , 47 r/ 3 ). Using the single coordinate, x = ki = 
^ 2 , we simplify (A 2 ) and (A5) to read 

5 (x) = l-be“-be*2^, (A 6 ) 

f{x) = 4 sin X — 2 sin 2x, (A7) 

and obtain the scaled band energies 


e± 


± 



2 


1/2 


f{x) + \g{x)\^ 


(AS) 


Here we define k = (fci/27r) bi + (^ 2 / 271 ) 621 with bi 2 The two bands are mirror-symmetric with respect to the 
denoting the reciprocal lattice vectors, and (fci,fc 2 ) G line e = 0 , and the dispersion of the upper band is shown 
























11 


in Fig. 7 for a few selected values of the governing param¬ 
eter F = JijJ\. Weak NNN hopping amplitudes open 
small topological gaps at the Dirac points, as illustrated 
by the dotted black line corresponding to F = 0.02. 
Stronger couplings produce relatively dispersionless por¬ 
tions of the bands. In Fig. 7, this result is illustrated by 
the full red and dashed blue lines drawn, respectively for 
Fq = l/-v/32 and Fj:, = l/\/27. The two reference values 
are obtained from different flatness criteria, and are very 
similar. The first value, Fq, is defined by the requirement 
that the Taylor expansion of the energy e+/J around the 
midpoint x = tt starts with the quartic (rather than the 
second order) term. The latter value, F/j, is the smallest 
at which the bandgap at the Dirac points becomes equal 
to the gap at the midpoint of the BZ diagonal. These 
two alternative criteria lead only to miniscule variations 
in the resulting manybody bandstructure, and we choose 
to adopt the value Fq as the reference. 


Appendix B: Operators in quasimomentum 
representation 


In the main text, we discuss the constituent parts of 
the effective Hamiltonian using their real-space forms 
which typically lend themselves to physically transpar¬ 
ent interpretation as combinations of hopping and inter¬ 
action events. The purpose of the present Appendix is 
to give the corresponding expressions in the reciprocal 
space, used in the actual numerical work. 

Denoting the number of elementary cells in the lattice 
As, we define the reciprocal-space creation operators 
and aj[,g in terms of the usual Fourier transforms 




^Ia 


t r 


2GA 


^at ifc (ri-6i) 

AT ® 




(Bl) 

(B2) 


iGB 


with ff(k) = ^i). The interaction Hamilto¬ 

nian reads for bosons 

r(b) 1 


Ant = ({fe}) al^A^LA^ksAak^A 


W 

Wb ({fe}) = ^ '^ki+k2,k3+k4’ 


(B4a) 

(B4b) 


and for fermions 

Ant = XI 4iB42A«fc3Aafc4B ■ Wf ({fe}), (B5a) 


{k} 

V 


({^}) — ^ 5(^1 ~ ^4) '^fci+fc2.fc3 + fc4’ 


(B5b) 


with the periodic Kronecker delta that allows the quasi¬ 
momenta in its arguments to differ by an elementary 
translation on the reciprocal lattice. The symbol {fe} = 
(fci, fc 2 , fea, A; 4 ) stands for the set of four quasimomenta 
involved in a scattering event. 

Proceeding to the second order we write 

A = -JJi(a) X [aj^B®fcAhBA(fc) + al^akBhAB(k), 

k 

(B6) 

with 


3 

/iBA(fe) = X (B7a) 

3 

/iAB(fe) = - X (B7b) 


and evaluating the commutator in Eq. (18b) obtain 



1 

hu! 



V3J"Ai"(a) 

fujj 


X X^(^) 






(B8) 


with the conjugate version of this equation applicable to 
annihilation operators. Note that in the quasimomentum 
representation we explicitly specify the sublattice index 
A or B while in the real space this was avoided by using 
the shorthand notation i versus f + /i. 

The first-order contribution to the effective Hamilto¬ 
nian is now written 


with f{k) = 2^ sin A; • Uj. These terms are diagonal 
and describe NNN hopping on separate sublattices. 

To describe the coupling of micromotion and interac¬ 
tions, we evaluate the nested commutators in Eq. (18c) 
and obtain, depending on the statistics, 

^(3.b|f) ^ _A^^^vpB|f({fe})5b|f, (B9) 


= Hint - JMa)' 


^ol^ttkAgik) + h.c. , (B3) 


with 
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‘5b({fe}) — + tlABiks) /iAB(^3) ®fciA®l2A®fc3AQfc4A ^ab(^i) ®IiB®fc2A®^3B®fc4A 

- hABi^s) h*AB{f«2) 4iA42B«fc3Bafc4A + hABi^s) h^Aif^i) ol^A^hA^k^Bak^B 

+ ^AB(fe4) ^Ab(^ 4) flfcjA®fe2A®^3Aflfc4A ~ ^abO^a) ^Ab(^i) ®feiB®l2A®^3A®fe4B 
~ ^Ab(^4) ^Ab(^ 2) ®IiA®l2B®fc3AOfe4B + ^Ab(^4) ^Ba(^ 3) QfcjA®l2A®fe3BQfc4B 

+ /iBA(fci) h*ABi.k2) 4iBaLB«fe3Aafc4A “ /iBA(fci) h^Aika) ai^B^i^A^ksBak^A (BIO) 

~ ^BA(fel) ^Ba(^ 4) ®fciB®l2A®fc3AOfe4B + ^BA(fcl) ^Ba(^i) A®l2A®fe3Aflfc4A 
+ ^BA(fc2) ^Ab(^i) ®fciB®l2B®fe3AS.fc4A ~ ^BA(fc2) ^Ba(^ 3) A®l2B®fc3B®fe4A 
~ ^Ba(^2) ^Ba(^ 4) ®fc4A®l2B®^3ASfe4B + hBA{k2) ^Ba(^ 2) flfc4A®l2A®fe3Aflfc4Ai 
“t” /z.C. “t" ‘ 1 ^'^ 

‘5f({fc}) = - /lBA(fe4) hXB{k2) ai^B^i^B^kiAakiA - hBA{k2) h*AB{k4) ai^B^i^B^kiAak^A 
~ kBAi^A) ^Ba(^i) ®IiA®l2A®fc3Aflfc4A + ^BA(fc4) ^Ba(^ 3) ®l4B®l2A®fe3BOfe4A 
+ ^BA(fe4) ^Ba(^ 4) al4B®l2A®fc3Aafc4B + ^BA(fe2) ^Ba(^i) ®fciA“fc2B®fc3Aafc4B 
+ ^BA(fe2) ^Ba(^ 2) ®l4B®l2A®^3ASfc4B ~ ^BA(fc2) ^Ba(^ 3) ®fc4B®fc2B®fe3B®fe4B 

- /lAB(fe3) h*AB{k2) ai^B^i^B^kiBak^B + hABika) h^Biks) ai^B^i^A^ksACtkiB (Bll) 

+ hABiki) /iab(^ 4) al4B®l2A®fc3Bafc4A + ^AB(fel) ^Ab(^i) ®fciB®fc2A®fc3Aafe4B 
+ ^AB(fel) ^Ab(^ 2) ®l4A®l2B®fc3A®fc4B ~ ^AB(fcl) ^Ab(^ 4) Si[.^^aj[^^afc3ASfc4A 

- hABiks) /iBA(fcl) aIiA®LA®fc3Bafc4B - /lAB(fcl) ^BAiks) al^A^k2A^k3BakiB 

+ h.c.. 
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